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^ The primordial abundances of light elements produced in the standard theory of Big Bang 
^ . nucleosynthesis (BBN) depend only on the cosmic ratio of baryons to photons, a quantity 
inferred from observations of the microwave background.^ The predicted^ primordial ^Li 
abundance is four times that measured in the atmospheres of Galactic halo stars.'^This dis- 
crepancy could be caused by modification of surface lithium abundances during the stars' 
'pj lifetimes^ or by physics beyond the Standard Model that impacts early nucleosynthesisJ^^iSl 
The lithium abundance of low-metallicity gas provides an alternative constraint on the pri- 
mordial abundance and cosmic evolution of lithiumPl that is not susceptible to the in situ 
modifications that may affect stellar atmospheres. Here we report observations of interstel- 
lar ^Li in the low-metallicity gas of the Small Magellanic Cloud, a nearby galaxy with one 
quarter of the Sun's metallicity. The present-day ^Li abundance in the Small Magellanic 
Cloud is nearly equal to the BBN predictions, severely constraining the amount of possible 
subsequent enrichment of the gas by stellar and cosmic ray nucleosynthesis. Our measure- 
ments can be reconciled with standard BBN with an extremely fine-tuned depletion of stellar 
Li with metallicity. They are also consistent with non-standard BBN. 

We obtained high resolution spectra (R ^ 70, 000) of the star Sk 143 (AzV 456), an 09.5 
lb star in the Small Magellanic Cloud (SMC), using UVES^^ on the 8.2-m VLT (observational 
details are given in the Supplementary Information). The sight line to this star was chosen for 
observation because it shows significant absorption from neutral atoms and molecule^^^'^ and a 
weak interstellar radiation field,'^ all of which favor the presence of neutral lithium, Li I. Li I 
absorption is clearly detected along this sight line (Figure HI)- 
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The derivation of the total Li/H abundance in the interstellar medium (ISM) requires large 
corrections for ionization, since A^(Li) ^ A^(Li II) ^ A^(Li I), and for the incorporation of Li into 
interstellar dust grainsplOur first approach to these corrections uses observations of adjacent ion- 
ization states of other metals, in this case Ca and Fe, to estimate the amount of unseen gas phase 
lithium. Assuming ionization balance and only atomic processes, the ratio N(Li Il)/A^(Li I) oc 
A^(Ca II) /N{Ca I) or oc A^(Fe II) /A^(Fe I), where the constant of proportionality involves the ra- 
tios of ionization rates and recombination coefficients for the elements in questionpHSl xhe ratio 
of ^Li I to total hydrogen in the SMC is log[A^(^Li l)/iV(H)] = -11.17 ± 0.04 (all uncertainties 
are la unless noted), where A^(H) = A^(H I) + 2A^(H2). Applying ionization corrections derived 
from Ca and Fe yields logarithmic abundances A^U) = log[iV(^Li)/iV(H)] + 12 = 2.79 ± 0.11 
and 3.01 ± 0.12. These calculations do not include more complicated (and uncertain) effects such 
as grain- assisted recombination,!^^!!!! nor do they correct for dust depletion. 

Our second approach uses the observatiorP' that A^(^Li l) /A^(K I) along sight lines through 
the Milky Way is nearly constant (with new determinations giving consistent result^^'^. When 
a differential ionization correction is applied, ''Li/K in the Milky Way ISM is consistent with the 
solar system ratio. Thus, ^Li and K appear to have very similar ionization and dust depletion be- 
haviors, and ^Li I/K I gives a good measure of the total (gas-i-dust phase) ^Li/ K.™!!] We measure 
log[A^(^Li I) /N{K I)] = -2.27±0.03 in the SMC, in agreement with the Galactic relationship .1™ 
Applying an ionization correction of +0.54 ± 0.08 dexP^ gives log[A^C^Li)/A^(K)] = -1.78 ± 
0.09. With the solar system ratio log(^Li/K)0 = —1.82 ± 0.05 derived from meteorites]^ we 
find [^Li/K]sMC = log[A^(^Li)/iV(K)] - log(^Li/K)Q = +0.04 ± 0.10. The ratio of ^Li to metal 
nuclei in the SMC is consistent with that found in the solar system and the Milky Way ISM^^: 
(7Li/K)sMC ~ ('Li/K)©. 

Although the ionization and depletion characteristics of S I are not as well tied to those of 
Li 1,1^ a similar approach using S I yields [^Li/S]sMC = —0.26 ± 0.11. The sub-solar ratio is 
consistent with a modest (0.3 dex) depletion of Li and K onto dust in the ISM^^ relative to S. 

We estimate A(^Li) by scaling ^Li/K to Li/H: A(^Li)sMC = ^('Li)0+[Fe/H]gMc+[K/Fe]gMc+ 
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pLi/K]sMc- We adopt fLi/K]sMc from above, the meteoritic A(^U)q = 3.23 ± O.OSPwith a 
mean present-day SMC metallicity [Fe/H]sMC = -0.59 ± 0.06 and an SMC K/Fe abundance 
[K/Fe]sMC = +0.00 ±0.10 (the last two discussed in the Supplementary Information). This yields 
A(^Li)sMC = 2.68 ± 0.16. Similarly scaling the ^Li/S resuk gives 2.38 ± 0.17. 

Most previous observational constraints on the primordial Li abundance have relied on mea- 
surements of atmospheric abundances in low-metallicity Galactic stars. Our first detection of in- 
terstellar lithium beyond the Milky Way opens a new window on the lithium problem. While 
there are significant uncertainties associated with ionization and dust effects, as demonstrated by 
the significant spread in A(^Li)sMc values, these are largely independent of the uncertainties that 
might affect stellar measurements of the primordial lithium abundance. Our recommended abso- 
lute abundance is A(^Li)sMc = 2.68 ± 0.16, or (^Li/H)sMC = (4.8 ± 1.8) x 10"^°, derived from 
''Li/K. This is compared to stellar ^Li abundanced^E^ at different metallicities in Figure [21 Our 
best estimate overlaps the prediction from standard BBN using the baryonic density deduced from 
the 5-year WMAP data^ ^(^Li) = 2.72 ± 0.06 (95% c.l.)^ although this leaves little room for 
the post-BBN chemical evolution,!^^*^ i.e., the contribution of freshly-synthesized Li to the ISM 
by stellar and cosmic ray nucleosynthesis (see representative model^in Figure [21). Our estimate 
of A(^Li)sMC is also consistent with the upper envelope of Li abundances in Milky Way thin disc 
stars (Figure [2l).ISl 

However, given the uncertainties in scaling to Ai^hVjsMC, the stronger result is our measure- 
ment [^Li/K]sMC = +0.04 ± 0.10. We compare [^Li/K]sMc with measurement^™ of fLi/Fe] 
and chemical evolution model^in Figure [3l The stars show a rapid decrease in [''Li/Fe] with in- 
creasing metallicity until [Fe/H] ^ —1, at which point the Li abundance increases roughly in lock- 
step with Fe such that disc stars have a nearly-constant [ ''Li/Fe] ratio similar to that found in the so- 
lar system. Our measurement of the present day ^Li-to-metal ratio in the SMC is in agreement with 
the nearly constant values found in the atmospheres of Milky Way disc stars (— 1 < [Fe/H] < 0), 
most of which formed >4 Gyr ago, with the solar system, and the modern-day Milky Way ISM.I^ 

Both the thin disc stars and our SMC measurements are below standard BBN predictions with 
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reasonable assumptions about post-BBN production, although it is often assumed these stars have 
had significant depletion of their surface Li abundance.!^ Taken at face value, the consistency of 
our SMC measurement with the [^Li/Fe] for those stars calls this assumption into question. While 
the models in Figures |2] and [3] are imprecise given the uncertain Li yields from stellar sources, 
they illustrate the tension between standard BBN predictions and our measurements if there is 
any post-BBN Li production. This tension can be relieved if a metallicity-dependent depletion 
of Li in stellar atmospheres is fine tuned in such a way that it is very strong below [Fe/H] ^ 
[Fe/H]sMC = —0.6 (to create the Spite plateau and avoid overproducing Li in the SMC ISM) and 
negligible at or above the SMC metallicity, conspiring to create a constant [''Li/Fe] ratio above 
[Fe/H] ^ —1. Alternatively, non-standard BBN scenarios can be invoked to allow for a lower 
primordial Li abundance.'^'^ 

If non-standard Li production occurs in the BBN epoch, many such models predict excess 
^Li compared with standard BBN. The only known source of post-Big Bang ^Li is production via 
cosmic ray interactions with ISM particles. Excess ^Li at the metallicity of the SMC would support 
non-standard production mechanisms, either in the BBN epochP^lor through the interaction of pre- 
galactic cosmic rays with intergalactic helium. Measurements of ^Li in stellar atmospheres are 
extremely difficult since the stellar line broadening is well in excess of the isotope shift. However, 
the ^Li I doublet is well-separated in our data due to the very low broadening in the cool ISM 
probed by Li I absorption. Our best fit to the SMC Li I absorption gives (^Li/''Li)sMC = 0.13 ± 
0.05 (see Supplementary Information and Figure [T]), giving a formal limit to the isotopic ratio in 
the SMC of (^Li/''Li)sMc < 0.28 (3cr). With higher signal-to-noise and resolution it should be 
possible to lower the limits for the interstellar isotope ratio in the SMC to a point that they provide 
constraints on non-standard BBN models. This approach has the advantage that ionization and 
dust depletion effects are not important for comparing the two isotopes of Li,^' making ^Li/^Li a 
powerful diagnostic of nucleosynthesis and non-standard evolution of Li abundances. 
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Figure 1 : Interstellar absorption by several neutral species seen toward Sk 143. Normalized interstellar 
absorption profiles from UVES plotted versus Local Standard of Rest (LSR) velocity and profile fit of the 
Li I absorption. The empirically-determined SNR is ~ 275 per pixel (5 pixels per resolution element) for 
the Li I observations. The full set of optical and ultraviolet absorption profiles seen toward this star and 
column densities measured from these are given in the Supplementary Information. The lower three panels 
show the profiles of Li I, K L and Fe I; the SMC cloud bearing Li I at tjlsr ~ +121 km is marked with 
the dashed line. The grey regions near Li I are possibly contaminated by diffuse interstellar bands or residual 
fringing, which may extend into the region containing Li absorption. The effects on the ^Li I columns are 
within the quoted uncertainties. The Li I absorption is composed of (hyper)fine structure components of both 
''Li I and ^Li I (shown respectively by the green and blue ticks in the top panel). The strong line along of 
''Li I is detected with ~ 16cr significance in the ISM of the SMC. A model fit to the Li I absorption complex 
is shown in the top panel (see Supplementary Information), with the fit residuals shown immediately below 
(normalized to the local error array). The free parameters for the fit are the polynomial coefficients for the 
stellar continuum, the central velocity, Doppler parameter (6-value), and column densities of ^Li I and ^Li I 
for the interstellar cloud. The red curve shows the best fit joint fit including both ^Li I and ^Li I, shown 
in green and blue, respectively. The best fit isotopic ratio is A^(^Li I)/A^(''Li I) = 0.13 ± 0.05 (68% c.l.), 
consistent with the presence of ^Li along the sight line, although below the Scr detection threshold. 
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Figure 2: Estimates of the lithium abundance in the SMC interstellar medium and in several different 
environments. Our best estimate for interstellar gas+dust phase abundance of ^(^Li) in the SMC is shown 
as the red circle with black core derived from the ^Li I/K I ratio. The present day metallicity of the SMC 
from early-type stars is [Fe/H] = —0.59 ± 0.06. (All uncertainties are Icr.) The point marked BBN 
and dotted horizontal line show the primordial abundance predicted by standard The green curves 

show recent model^ for post-BBN ^Li nucleosynthesis due to cosmic rays (CRs) and stars. By adjusting 
the yields from low-mass stars, the models are forced to match the solar system meteoritic abundance'^ (see 
Supplementary Information). The solid and dashed lines correspond to models A and B^^ which respectively 
include or not a presumed contribution to ^Li from core-collapse supernovae. The blue hatched area shows 
the range of abundances derived for Population II stars in the Galactic haloP with the "Spite plateau" in 
this sample at A(^Li)popii ~ 2.10 ± O.lO.'^The violet hatched region shows the range of measurements 
seen in Galactic thin disk stars, where the thicker lines denote the six most Li-rich stars in a series of eight 
metallicity bins.^The selection of thin disk stars includes objects over a range of masses and temperatures, 
including stars that are expected to have destroyed a fair fraction of their Li. Thus, the upper envelope of 
the distribution represents the best estimate of the intrinsic ISM Li abundance at the epoch of formation for 
those stars, and the thicker dashed lines for the thin disk sample are most appropriate for comparison with 
the SMC value. The most Li-rich stars in the Milky Way thin discP^ within 0. 1 dex of the SMC metaUicity 
give A(^'Li)mw = 2.54 it 0.05, consistent with our estimate A(^Li)sMc = 2.68 ± 0.16. 
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Figure 3: Estimates of Li/Fe in the SMC interstellar medium and in several different environ- 
ments. The SMC value is derived from the ^Li I/K I ratio. At low metaUicities ([Fe/H] < —1), stellar 
measurements^ trace the build-up of Fe with a constant Li abundance along the Spite plateau. At higher 
metaUicities, disc star abundances^^ show a turnover to roughly constant [^Li/Fe] at values consistent with 
the solar system/meteoritic valuPl (shown as the dash-dotted Une). Our SMC estimate is consistent with 
the solar system and disc star abundances in this region of relatively constant ''Li/Fe abundances, with 
[^Li/Fe]sMC = +0.04 ± 0.14 for [K/Fe]sMC = 0.0 ± 0.10 (Supplemental Information). The most Li-rich 
disc stars within 0. 1 dex of the SMC metallicity have a mean Li/Fe]) = -0.13 ±0.05. (All uncertainties 
are la.) The green curves show the chemical evolution model^as in Figure[2l while the dotted line shows 
the behavior of [^Li/Fe] for the standard BBN primordial abundance with no subsequent evolution of ^Li. 
The relative uniformity of the stellar ^Li/Fe abundances at [Fe/H] > — 1 could be caused by a delicate 
balance of Li and Fe production and metallicity-dependent Li astration (not ruled out given the changes 
in mean age and mass potentially present in the sampl^^^. However, the agreement of the [^Li/Fe] ratio 
seen in these old stars (ages >4 Gyi^^ and in the present-day interstellar medium of the SMC suggests 
little change in the stellar abundances for metaUicities [Fe/H] > —0.6 through the solar metallicity. To 
bring the stellar and SMC interstellar abundances into agreement with standard BBN predictions requires 
a delayed injection of significant ^Li from stellar production mechanisms as well as vigorous depletion of 
stellar surface ^Li abundances at metaUicities just below that of the SMC. 
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Supplementary Information 



1 Observations and Data Reduction 

We use spectroscopic data from two instruments in our analysis: ground-based optical observations 
from the Ultraviolet Echelle Spectrograph (UVESj^on ESO's Very Large Telescope (UT-2), and 
space-based ultraviolet observations from the Space Telescope Imaging Spectrograph (STIS) on- 
board the Hubble Space Telescope. 

UVES Data The UVES data presented here were taken in service mode on 2008 September 25 
under the program 382.B-0556(A). The observations consist of 10 exposures of 2850 seconds each. 
We used the 0'.'7 slit, feeding the light to dichroic #2. The grating central wavelengths were 3,900 
A in the blue and 7,600 A in the red. No binning was done on the chips. The data were taken 
through thin clouds and with seeing ranging from ^ I'.'O to I'.'S. 

The data were reduced using the UVES pipeline provided by ESO. In order to achieve high 
signal-to-noise, we used the pixel-by-pixel flat field correction, setting "FFMETHOD=pixel" in 
the reduction pipeline and using optimal extraction of the data. The individual exposures were 
corrected for heliocentric motion and coadded, weighted by the variance. The signal-to-noise ratio 
(SNR) varies from SNR ~ 150 per pixel in the far blue (e.g., near Ti II 3,383) to ~ 300 in some 
regions of the spectrum (e.g., the Ca I and CH+ transitions near 4,200 A). Near the Li I transition 
at A ^ 6, 710 A we measure SNR ~ 275. 

UVES Line Spread Function: The width of the UVES line spread function (LSF) is important for 
deriving accurate fit parameters for the blended isotopes of Li. While ultimately the column density 
of Li I is not strongly changed by small changes in the LSF width, the 6-values are somewhat 
sensitive to the LSF width. The UVES pipeline produces estimates of the width of the LSF on 
the basis of observations of ThAr lamps taken for wavelength calibration. For the blue CCD data, 
the ESO quality-control summaries available on-line suggest the resolution for our set-up during 
this time is i? ^ 66,000 or At> ^ 4.54 kms~^, which is appropriate for roughly the month 
surrounding our observations. The ThAr lines in our own calibration observations yield a value 
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indistinguishable from this. 

On the red chip (REDL), notably in the vicinity of the Li I lines near 6,710 A, the mean 
FWHM of these lines from our calibration data is 4.47 km s^^ for the 38 lines in the range 
6, 700 A < A < 6, 720 A and 4.42 km s"^ for the 225 lines over the broader range 6, 650 A < A < 
6, 770 A. These correspond to a resolution R ^ 68, 000. 

However, in the course of fitting the Li I blend, adopting this resolution seemed to under- 
predict the cores of the lines and provide too much absorption immediately outside of the cores, 
suggesting the adopted LSF may be too broad. D. Welty kindly provided us with Na I 5,889/5,895 
observations of the Sk 143 sight line taken with the CES spectrograph on the ESO 3.6-m telescope 
with a resolution Av = 1.35 km (FWHM) or R ^ 222, 000. These data were presented in 
Welty et al.'^ We derived the best LSF breadth for the current UVES red CCD data by compar- 
ing smoothed CES data with the UVES observations of these absorption lines, using the as an 
estimate of the goodness of fit. The free parameters of the model in this case were the breadth of 
the smoothing kernel (assumed Gaussian) and the background offset of the CES observations (for 
which the original background subtraction was derived by matching equivalent widths of the Na I 
lines with earlier UVES data-Sl). 

The final result of this approach yields a best fit LSF breadth of Av = 4.03 ± 0.08 km 
(FWHM) or resolution R = 74, 300, which we adopt in fitting the Li I blend. The difference 
between our derived LSF and the values from the ThAr emission lines is likely due to the manner 
in which the light fills (or not) the slit during the observations. We caution that we assume in our 
fitting a single Gaussian LSF; the true LSF may be more complex than this. However, we have 
tested the impact of the LSF on our measurement of both the ''Li I column and isotopic ratio (see 
below), and find the results robust to even 5a variations in the LSF breadth. This LSF breadth likely 
only applies to the data taken with the REDL chip. LSF changes with wavelength and detector are 
to be expected. However, the breadth of the LSF for the data taken with the blue side may be 
smaller than the 4.55 km s^^ predicted from the ThAr exposures. 
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STIS Data We make use of archival STIS observations taken with the E140H and E230H gratings 
using the 0.20" x 0.09" apertures. These data were acquired under program 9383 (PI: K. Gordon) 
and have been previously reported in the literature. The resolution of these data is i? ~ 114, 000 
corresponding to Av ^ 2.75 km s^^. The data were processed with CALSTIS v2.27 and otherwise 
reduced following the discussion in earlier papers.l^ The SNR of these data are limited by the 
relatively low UV flux of this highly extincted star for high resolution spectroscopy and ranges 
between 5 and 15 (per pixel) for the lines of interest. 

2 Hydrogen column density along the Sk 143 sight line 

We derive the H I column density toward Sk 143 by fitting the strong Lya absorption seen in the 
Hubble Space Telescope (//57)/STIS E140H observations. Following Lehner et al.P^ we fitted Lya 
with two components (Galactic and SMC) at fixed velocities derived from the metal absorption 
lines and the H I 21 -cm emission profile for the Galactic component (see below), which is in 
agreement with the velocity from the absorption profiles. Figure |ST] shows the STIS data for the 
region surrounding Lya (binned by 10 pixels for display purposes) along with our best fit H I 
profile and continuum. For the SMC component we find log A^(H I) = 21.07 ± 0.05. We tested 
the robustness of our fit by undertaking several simulations wherein the velocity centroids of the 
gas were changed by ±5 km s^^, the placement of the continuum was varied, and the continuum 
was modeled by a range of polynomials of degree of 2 to 4. The range of values derived in these 
simulations is used as an estimate of the errors associated with these effects and is included in 
our final uncertainty. The column density of the Galactic component in the fit is quite uncertain 
because it is much weaker (see 21-cm emission spectrum discussed below and in Figure [S2l). Its 
weakness is in fact a benefit for deriving an accurate SMC A^(H I) as the damping wings are 
dominated by the SMC component. 

Table [ST] summarizes previous determinations of the H I column density in the SMC com- 
ponent toward Sk 143. Previous estimates relied on H I 21-cm emission data or Lya absorption 
from a low resolution and low signal-to-noise spectrum from the International Ultraviolet Explorer 
{IUE)^^^T\vQ poor quality of the available lUE spectrum would seem to preclude a precise mea- 
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Figure SI: Fit to interstellar H I Lya absorption observed by HST/STIS against the stellar continuum 

of Sk 143. The hlue line shows the adopted continuum, while the red line shows the continuum absorbed by 
an H 1 column density of log A^(H 1) = 21.07 ± 0.05. 
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Figure S2: H I 21-cm emission line observations of the sight line to Sk 143. The emission line data in 
black are taken from the Leiden/Argentine/Bonn (LAB) Survey'^, with a 36' beam, while the data in red are 
taken from the Parkes/ATCA survey of the SMC——, with a 1 .6' beam. The principal SMC component as 
seen in absorption is marked. 
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surement, and the first estimate using those data had a large uncertainty of 40-50%P"^ Most 
recent studies of this sight linePE^H have adopted H I estimates derived from the /f/£'-observed 
Lya absorptiorP"^, or analyses of the H I 21 -cm emission line observations.'^'^ The most recent 
determinations using lUE^^ give high values of log A^(H I) ^ 21.6, consistent with values derived 
by integrating the full SMC profile of the H I 21 -cm emission seen in Figure [S2l For comparison 
with absorption line data toward Sk 143, however, the H I 21 -cm emission estimate is not always 
ideal, as it includes contributions from gas behind the star. Most of the absorption toward Sk 143 
is found at vlsr ~ 121 km s^^ as seen in Figures [S3] and [SH This component is seen as well in 
21-cm emission, but the 21-cm profile contains significant amounts of gas to much larger veloc- 
ities. Some of this is seen in the absorption line observations, but with a smaller contribution to 
the total column of undepleted elements (e.g., see the total gas column traced by S II absorption). 
Integrating the combined Parkes and ATCA H I profileP^ in Figure [S2] over the velocity range 
containing the majority of the optical/UV metal absorption, vlsr = 80 to 146 km s^^, we find 
logA^(H I) = 21.17 ± 0.11 (where the error is a combination of statistical and "beam" errors). 
Thus, the integration over this smaller velocity range gives results consistent with our measure- 
ment and another recent one of Welty & Crowthei'^^ at < 2a. We note, though, that even though 
there is a kinematic correspondence in the component at 121 km in the absorption and emis- 
sion profiles, there is a possible (unknown) uncertainty owing to the unknown depth of the star in 
the SMC. The moderate difference between our result and that of Welty & Crowther is due to the 
difference in the H I column assigned to the Milky Way, for which they derive values from 21-cm 
measurements while we fit the Lya directly. 

In what follows, we adopt our H I column density measurement from the Lya fitting. With 
the H2 column derived from FUSE data,'^' the total hydrogen column along this sight line is then 
log A^(H) = log[iV(H I) + 2A^(H2)] = 21.46 ± 0.04. The sight line has a high molecular fraction, 
with /(H2) = 2iV(H2)/iV(H) = 0.6. 
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Table S 1 : H I column density estimates to Sk 143 using Lya 



Reference 


logAf(HI)sMc 


Data 


This Work 


21.07 ±0.05 


HST/STIS 


Welty & Crowther (2010)^^ 


20.95 


HST/STIS 


Gordon et al. (lOOSp^ 


21.60 ±0.05 


WE 


Bouchet et al. (1985jSl 


21.40 ±0.17 


WE 


Fitzpatrick(1985jSl 


21.18±0.18 


WE 



3 Metal column densities along the Sk 143 sight line 

Table |S2] gives our adopted final column density for all of the atomic and ionic species detected to- 
ward Sk 143 along with the methodology used to derive the columns. We measure metal atom and 
ion column densities along the Sk 143 sight line using a combination of profile fitting and apparent 
optical depth integration methods, depending on the quality of the data and the specifics of the 
individual absorption lines. We assume oscillator strengths from the compilation of Morton. For 
each atomic or ionic species, two column densities are given: the first applies to the column den- 
sity integrated over the entirety of SMC absorption, the second gives the column density associated 
with the +121 km s^^ component within which the Li I and other neutral species arise. In some 
cases (e.g., Li I) these are the same as no significant absorption exists outside of this cold, dense 
component. We use the column densities for this principal component at wlsr ~ +121 km s^^ in 
all calculations in the main text, since the Li I-bearing gas is all contained in this component. 

Our component fitting column densities are derived using the software package 
made available by R. Carswell ('http://www.ast.cam.ac.uk/~rfc/vpfit.html). VPFIT minimizes the 
goodness of fit parameter using a Gauss-Newton type system for parameter updates."^" The free 
parameters of the fit for each assumed interstellar "component" are the central velocity, Doppler 
parameter (6- value), and column density of the component. We fit each species independently, with 
no assumptions about 6-values or central velocities from other absorbing species (as opposed to 
assuming the 6-values are linked through a common temperature and non-thermal velocity). This 



6 



Table S2: Column densities toward Sk 143 



Species 


log N{SMC) 


log7V(+121 km s-i) 


Method 




[X/H]sMC 


HI 


21.07 ±0.05 


21.07 ±0.05 


1 


= 12.00 




H2 


20.93 ± 0.09 


20.93 ±0.09 


2 


= 12.00 




Htotal 


21.46 ±0.04 


21.46 ±0.04 




= 12.00 




^Lil 


10.29 ±0.02 


10.29 ±0.02 


3 


3.23 ±0.05 


-0.55 ±0.16 


^Lil 


9.41 ±0.15 


9.41 ±0.15 


3 


2.14 ±0.05 


-0.38 ±0.21 


Nal 


12.96 ±0.04 


12.96 ±0.04 


3 


6.26 ±0.02 


> -2.76 ±0.06 


Mgl 


> 13.11 ±0.06 


> 13.11 ±0.06 


4 


7.53 ±0.01 




Mgll 


15.37 ±0.09 


15.37 ±0.09 


3 


7.53 ±0.01 


-1.62 ±0.10 


SI 


14.15 ±0.08 


14.15 ±0.08 


3 


7.14 ±0.02 




S II 


> 15.36 ±0.06 


> 15.17 ±0.05 


3 


7.14 ±0.02 


> -1.24 ±0.07 


KI 


12.61 ±0.02 


12.61 ±0.02 


3 


5.08 ± 0.02 


> -1.93 ±0.04 


Cal 


9.42 ±0.10 


9.42 ±0.10 


4 


6.30 ±0.02 




Call 


12.42 ±0.01 


11.68 ±0.01 


3 


6.30 ±0.02 


> -3.34 ±0.04 


Fel 


11.83 ±0.05 


11.83 ±0.05 


4 


7.45 ±0.01 




Fell 


14.63 ±0.07 


14.15 ±0.10 


4 


7.45 ±0.01 


-2.28 ±0.08 


Till 


11.88 ±0.03 


11.47 ±0.03 


4 


4.92 ± 0.03 


-2.50 ±0.05 


Ni II 


13.101H? 


.2 07+0.08 


4 


6.20 ±0.01 


-2.56 ±0.11 


Zn II 


> 12.80 ±0.06 


> 12.80 ±0.06 


4 


4.60 ± 0.03 


> -1.26 ±0.07 



Methods used for determining column densities: (1): Lya Fitting; (2): Fitting by Cartledge et al.pl(3): 
Profile fitting; (4): Na{v) integration. The values of [X/H]gMC are gas-phase abundances, i.e., they reflect 
the sub-solar SMC metalhcity and the differential incorporation of elements into the sohd phase (dust 
grains). 
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has the advantage that additional unsuspected components should not adversely affect the results, 
so long as the lines are not strongly saturated. We fit multiplets simultaneously where available. 
We did not, however, simultaneously fit the strongest doublets of Na I and K I with the weakest 
doublets available (see Figures [S3] and [S4l for the observed lines). The column densities in Table 
[S2]are derived based on the weaker doublets of these species due to the extreme saturation of the 
strongest lines. 

Our fitting of all of the neutral species suggests that the SMC absorption seen in weak lines of 
neutral species (i.e., not including the Na I D lines) are well described by a single component. We 
see no compelling evidence for additional components that contribute significantly to the column 
densities of the neutral species. The high resolution CES data presented by Welty et al.^^ are 
consistent with a single dominant component with only very weak additional components at higher 
relative velocities. We note that the presence of additional unknown components could affect the 
interpretation of the Li isotopic ratio if they had strong variations in the ratio between components 
or if there were moderate strength components present with separations matching the isotope shift. 
For derivation of total column densities of ^Li (and largely even ^Li), the presence of multiple 
components is not particularly important. In general our profile fits to these neutral species gave 
6-values of order 0.7 km consistent with previous studies of this sight line.'^^'^ We note also 
that the presence of additional unidentified components is unlikely unless they have a separation 
less than this 0.7 km s^^. 

Our derived column densities of Na I, K I, and Ca II vary somewhat from these earlier 
workd^^^ (with our columns being higher) because we have generally studied weaker lines, fully 
fit multiple lines, or in the case of Ca II used a careful treatment of the continuum for the weak line, 
which lies on the stellar He absorption requiring special techniques for the continuum assessment. 

Our fitting assumed Gaussian LSFs with breadths derived as described above for the UVES 
observations of Li I. For the blue CCD data we had no higher resolution comparison observations, 
so we adopted the LSF breadth suggested by the calibration lamp data. There was no obvious 
indication that this LSF is too broad. For fitting the STIS E140H observations of Mg II, we 
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Figure S3: Normalized interstellar absorption profiles for optical metal line transitions toward Sk 143 
plotted versus LSR velocity. These data are from the VLT+UVES observations of Sk 143. Contaminating 
absorption from other species are marked in light gray. 
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Figure S4: Normalized interstellar absorption profiles for ultraviolet metal fine transitions toward Sk 
143 plotted versus LSR velocity. These data are from the HST/STIS E140H and E230H observations of Sk 
143. Contaminating absorption from other species are marked in light gray. Optical observations from the 
previous figure are shown at the top for comparison. 
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assumed a breadth of 2.75 km (FWHM). The STIS LSF is not strictly a Gaussian. In our case 
we were only interested in separating the SMC Mg II absorption from the contaminating Milky 
Way absorption. For these purposes and with the strength of the lines, the assumption of a Gaussian 
LSF will not adversely affect the results. 

We have used profile fitting to derive column densities for neutral species with small b- 
values in case unresolved saturation is present or where blending is a concern. For most of the 
ionic species we adopt column densities derived through an integration of the apparent column 
density, Na{v), profiles,'^ largely following methodologies in our earlier works. Saturation can be 
identified by comparing the Na{v) profiles from two lines of an absorbing species with significantly 
different oscillator strengths (/-values). The entries in Table [S2] listed as lower limits are those for 
which we have evidence for unresolved saturation. 

Table [S2] also gives the solar system abundanced^ for each element and the relative gas- 
phase abundance of the element in the ISM of the SMC, [X/H]smc- The latter quantity reflects 
the base subsolar metallicity of the present-day SMC, [Fe/H]sMC = —0.59 ± 0.06 derived from a 
straight mean of observations of early type stars,!^^^ and the effects of differential dust depletion 
of the elements.S3 Where two ionization states are measured, we only give the gas-phase values 
for the most abundant ion. Inequalities either reflect the non-dominant nature of the species (e.g., 
for the neutrals) or the probable presence of unresolved saturation. The abundances listed for the 
Li isotopes are drawn from the main text (i.e., adopting the abundances derived from scaling the 
Li I/K I ratio). 

Assessing the effects of dust depletion on the gas along this sight line is somewhat compli- 
cated by the lack of a quality measurement of an undepleted metal species. The [X/H] values in 
Table[S2]measure the effects of the sub-solar metallicity of the SMC and its non-solar relative abun- 
dance pattern in addition to dust depletion effects. Welty & Crowthei^^have recently estimated the 
intrinsic Ti abundance of the SMC based on early-type stars, deriving an estimate for the present 
day abundance of y4(Ti)sMc ~ 4.28. Our measurements of the gas-phase Ti II then suggest Ti is 
depleted by —1.86 dex relative to the intrinsic SMC abundance; sight lines in the Milky Way have 
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depletions —2 to —3 for similar E{B — V) or CN columns, but it is not dissimilar from sight lines 
in the Galactic disk with modest depletions. Defining the reference abundance can be difficult for 
assessing the depletion of an element. Another approach to understanding the level of elemental 
incorporation into grains is to measure the differential depletion between two metals. The least de- 
pleted species with a good measure of the total column density in the table is that of Mg, for which 
we have a good measure of the dominant ion Mg II. Comparing strongly-depleted metals to Mg 
in the SMC gas toward Sk 143 gives [Fe/Mg]sMc = -0.66 ± 0.13, [Ni/Mg]sMc = -0.94 ± 0.15, 
and [Ti/Mg]sMC = -0.88 ± 0.11. The uncertainties are large given the low SNR of the STIS 
data for deriving abundances, so comparing these data to Milky Way sight lines does not lead to 
great insight. These are consistent with relative gas-phase abundances seen in Milky Way disk gas 
for sight lines with low to modest depletion (at least for the disk) .^^Ell Measurements of Li in the 
Milky Way typically probe sight lines or clouds with stronger depletions than these values; the Li 
depletion effects should not as strong in this SMC gas as in the gas typically probed in the Milky 
Way. 

4 Lithium abundance estimates 

Several estimates for the gas-phase ^Li abundance are given in the text, which are summarized 
in Table [S3l Estimating the gas-phase ^Li abundance from measurements of A^(^Li l)/A^(H) re- 
quires the application of an ionization correction to account for the unseen ^Li, which is mostly in 
its singly ionized form in the ISM. Assuming ionization rate balance and only atomic processes, 
N{U) /N{U I) ^ A^(Li II) /A^(Li I) = r(Li I) [^^^(Li II, T)]-\ where r(Li I) is the photoioniza- 
tion rate and a (Li II, T) the recombination coefficient. (We generically refer to Li I here, which 
applies to either isotope or their sum.) At first glance it would seem that knowledge of F and Ue are 
crucial to estimating the total abundance. However, similar equations can be written for observed 
adjacent ionization states, notably for A^(Ca Il)/A^(Ca I) and A^(Fe Il)/A^(Fe I). Substituting 
derived from these ratios into the expression for Li gives, for example, 

iV(Li II) _ A^(Fe II) r(Li I) /r(Fe I) 
A^(Li I) ~ iV(Fe I) a(Li II, T) /a(Fe II, T) ' 
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Table S3: Summary of SMC lithium abundances 



Quantity 


Value 


Methodology 


^('^Ll)sMC 


2.68 ± 0.16 


Scaled from 'Li I/K I 


^(^Ll)sMC 


2.38 ± 0.17 


Scaled from 'Li I/S I 


^(^Li)sMC 


3.01 ±0.12 


Ionization correction from Ca Il/Ca I 


^(^Li)sMC 


2.79 ±0.11 


Ionization correction from Fe II/Fe I 


[7Li/K]sMC 


±0.04 ±0.10 


Differential ionization correction 


rLi/s]sMc 


-0.26 ±0.11 


Differential ionization correction 



All absolute and relative abundances of SMC Li derived using sev- 
eral different methods are given. The recommended values are those 
derived by comparison with K, i.e., A(^Li)sMC — 2.68 ± 0.16 and 
[^Li/K]sMC +0.04 ± 0.10. The former relies on the application of 
several scale factors, including a differential ionization correction, the 
metallicity of the SMC, and an estimated K/Fe abundance. The latter relies 
on a differential ionization correction, but one that is well founded in the 
observed Galactic relationship between Li I and K I. 

While the individual photoionization rates depend on the intensity of the interstellar radiation 
field, to first order the ratio of photoionization rates r(Li l)/r(Fe I) does not. What is crucial 
in this case is the shape of the radiation field over the range of energies encompassed by the 
ionization edges of these neutrals. The grossly similar temperature sensitivities for the recombi- 
nation coefficients mitigates the effect of this unknown, as well. The total abundance of Li is then 
log[A^(^Li)/A^(H)] = log[A^(^Li l)/A^(H)] - loga;(Li I), where the term - logx(Li I) is the ion- 
ization correction factor and x(Li I) = A^(Li l)/A^(Li) is the ionization fraction of Li in neutral 
form. 

As discussed in the main text, we derive ionization corrections for the Li/H comparison on 
the basis of the observed ratios iV(Ca Il)/A^(Ca I) and A^(Fe Il)/A^(Fe I). These give ionization 
corrections of — log a;(Li I) = ±1.96±0.10 and ±2.18±0.12 from the Ca and Fe ions, respectively, 
assuming standard relative ionization rates with specific interstellar radiation fields .'^'^ These are 
derived using column densities from the wlsr ~ ±121 km cloud only. The ratio of ^Li I to 
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total hydrogen in the SMC gas along this sight line is logfA^C^Li I) /iV(H)] = -11.17±0.04, giving 
log[A^(^Li)/A^(H)] = -9.21 ±0.11 and -8.99 ± 0.13. While we also have data covering adjacent 
ions of S and Mg, we have only limits for one ionization state of each of these elements that do not 
produce useful limits on the ionization corrections. 

The discrepancies between these calculations are in part due to errors in the calculated pho- 
to ionization rates and recombination coefficients. However, it is also likely that the neutral and 
singly-ionized Ca and Fe, which have different ionization potentials, photoionization cross sec- 
tions, and recombination coefficients, trace slightly different conditions within the cloud. The 
conditions traced by these species will also be different than those traced by Li I and Li II, limiting 
the precision of this approach. In fact, electron densities derived from various tracers can vary, 
with Ca I typically predicting significantly larger rif, than most.l^ In the SMC cloud probed here, 
however, the Ca I column seems to be unusually low, as indicated by the Ca I to K I ratio, for ex- 
ample, compared with Milky Way sight lines. The integrated A^(Ca II) /A^(Ca I) is higher than all 
measured values in the compilation of Welty et al.'^ (although some sight lines with nondetections 
of Ca I could have higher values), although the ratio in the wlsr ~ +121 km s^^ component is 
near the mean of those measurements. 

Using Li I/K I as a tracer of the total interstellar Li/K abundance follows a similar approach: 

N{U II) _ N{U I) r(Li I) /r(K I) 
A^(K II) ~ iV(K I) tt(Li II, T) /a(K II, T) ' 

In this case the ratio of relative ionization rates to recombination coefficients gives an ionization 
correction of +0.54 ± 0.08 (following SteigmaiP^ in adopting an error based on the range of re- 
combination coefficients derived from a wide range of temperatures). The Li I/K I ratio is powerful 
given the very nearly linear relationship between the column densities of these species in the Milky 
Way.ESEHWe assume the same physics applies in the SMC, and thus this ratio should be similarly 
useful there. This implies Li and K have very similar dust depletiorP"^ and ionization conditions 
that track each other well. Similarly, the K I/S I ratio is a power law with an index somewhat 
greater than unity.l^ The tight correlation with larger power law slope likely implies the two neu- 
trals trace similar physical conditions, but with different depletion characteristics. We make the 
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assumption that Li I also traces physical conditions similar to those of S I-bearing clouds. The 
ionization correction relating Li I/S I to Li/S is —0.32 ± 0.08 dex. We do not use the other neutral 
tracers for such a derivation of the Li to metal nuclei ratio because of the weaker correlations and 
evidence for differeng physical conditions between Li I and the remaining neutrals. 

Transforming [^Li/K]sMC to an estimate of A(^Li)sMC involves a number of scale factors. 
Our adopted mean metallicity of the SMC ([Fe/H]sMc) is described above. One other scale fac- 
tor is the abundance of K/Fe with respect to the solar system. The K/Fe ratio is not measured 
directly in the SMC. However, K has its origins in explosive oxygen burning in high mass stars 
and tends to behave nucleosynthetically like an a element.!^ Studies of early-type stars |^^*^^^ A 
supergiantsj^ and H II region^ suggest [S, Si/Fe] ^ in the present-day SMC. These studies 
can have significant uncertainties, and there is some dispersion in the results. We adopt a value 
[K/Fe]sMC = +0.0 ± 0.1 in our analysis (uncertainty from dispersion in the a to Fe dispersion in 
the SMC). We note that this is at odds with what is seen in the Milky Way. The mean of two sam- 
ples of Milky Way stars within ±0.1 dex of the SMC metallicity,™! [K/Fe]sMC = +0.27 ± 0.09 
(standard deviation). The discrepancy is due to the significantly different star formation and en- 
richment history of the SMC compared with the Milky Way.!^ 

Our ^Li abundance and ''Li/K ratio is significantly higher than that impliecP^ by the upper 
limit^§3*^or claimed detection^ toward SN1987A in the Large Magellanic Cloud (LMC). With 
an abundance roughly twice that of the SMC, the LMC should have a ^Li abundance at least as 
high as that of the SMC. The values reported for this sight line vary significantly, and we view them 
with significant skepticism in hindsight. Aside from the discrepancies in the Li I measurements, 
there are worries about saturation corrections for the K I reference. Furthermore, the sight line to 
SN1987A probes a significantly lower column density, more highly ionized portion of the ISM. 
This sight line does not favor the formation of atomic species such as Li I, making the ionization 
effects even more severe. Future observations of Li I absorption within the LMC should provide 
more clarity to the chemical evolution of Li in the Clouds. 
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5 Lithium Isotopic Ratio 

The Li isotope ratio can provide some constraints on the production of Li and any non-standard 
contributions to Li since ^Li is not expected to be produced in significant amounts in BBN.'" Our 
estimate of the ^Li/^Li relies on simultaneously fitting the hyperfine structure of the ^Li I and ^Li I 
absorption. We adopt wavelengths for the hyperfine levels from Sansonetti et and /-values 
summarized in Welty et al^ We assume the ISM absorption from the neutral species is well- 
characterized by a single component or cloud and that both isotopes have a common 6- value. The 
latter assumption is only strictly correct if the velocity dispersion is dominated by non-thermal 
motions (which is likely), although it has little practical effect given the weakness of the Li I 
absorption. We fit a stellar continuum to the data before the fitting process. The continuum in the 
region around Li I was well fit with a first-order Legendre polynomial. We assessed higher order 
fits using an F-test, but they did not provide a statistically significant improvement in the goodness- 
of-fit parameter over the region used to derive the continuum. In the profile fitting process, we 
allow the parameters of the continuum fit to vary during the minimization process. This allows 
VPFIT both to esimate the continuum fit objectively and to include an estimate of the uncertainties 
caused by continuum fitting into the final error budget. 

The column densities reported in Table [S2] are derived from these fits. The best fit Doppler 
parameter is 6 = 0.8 ± 0.5 km s^^. Because the absorption is so weak, the 5-value has a minor 
effect on the quality of the column density determination (as would multiple components). The 
fits were made with the UVES LSF as discussed above. The resolution we derive for the red side 
data are higher than typically advertised or adopted for the instrument. Our error budget includes 
a contribution derived by refitting the data with LSF breadths varied by ±la, but the resulting 
differences in the central value are minimal. Adopting a significantly broader FWHM for the LSF 
(e.g., Av ^ 4.45 km s^^ as suggested by the ThAr lines) has very little effect on the central 
value or uncertainties, although the goodness-of-fit parameter is larger. In this case, the profiles 
appear visibly broader than the data, with the peaks of the models underpredicting the observed 
absorption in the core and the wings overpredicting the observed absorption in the edges of the 
profiles. We have tested the impact of adding additional components, but VPFIT rejects them as 
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unnecessary as they do not significantly improve the fit. 

It is important to point out that the potential presence of ^Li I absorption is not simply a 
function of the manner in which we have proceeded with the profile fitting. The ^Li I contribution 
is hinted at from simple equivalent width measurements of the ''Li I doublet. Using our initial 
continuum (i.e., not the final one adopted above through the formal fitting process), the weak and 
strong lines of '^Li I give equivalent widths Wx = 3.77 ± 0.23 mA and 2.40 ± 0.22 mA. (We 
caution both the errors and equivalent widths derived here may not be equivalent to those reported 
through the profile fitting analysis due to the different manner in which the continuum placement 
uncertainties are included and the crude nature of the separation of the strong and weak lines in the 
direct integration of equivalent widths.) Assuming there is no saturation, a safe assumption given 
the weakness of the lines, these lines should have equivalent widths in a 2:1 ratio. The observed 
ratio of 1.57:1 is lower due to the blended contribution of ^Li I to the weak member of the ''Li I 
doublet. A simple estimate of the ^Li I column density can be derived by subtracting half of the 
''Li I strong line equivalent width from the weak line value. This yields log A^(^Li I) ^ 9.41 ±0.20, 
completely consistent with the value derived from the profile fitting. 

The derived isotopic ratio for the SMC gas toward Sk 143 is ^Li/'^'LisMC = 0.13 ± 0.05, 
which implies a limit ^Li/^LisMC < 0.28 (3cr). The limit itself does not constrain the ratio in an 
interesting way. It is, for example, consistent with the Solar System value, (^Li/^Li)© = 0.0787 ± 
0.0004,^^ and it is very similar to a number of measurements in Milky Way gas,'^'^ although 
there are clouds and sight lines with large variations from this value.l^ Kawanomoto et al.^^ have 
recently derived a mean for the Milky Way ISM of ('5Li/'^Li)Mw = 0.13 ± 0.04. The only known 
source of post-Big Bang ^Li is production via cosmic ray interactions with interstellar medium 
(ISM) particles, either through spallation or a + a fusion.l^^^For standard energy distributions of 
Galactic cosmic rays, these isotopes are produced in a ratio ^Li/''Li ^ 0.67, while non-standard 
assumptions can push this ratio to 0.5 or a bit lower.^^ For the standard production ratio, our central 
value implies that (19 ± 8)% of the ^Li in the ISM of the SMC has been produced via cosmic rays, 
with a limit of < 42% (3a). 
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Our data show Li I at < 3cr significance, and the reality of the Li I at the levels of our best 
fit value should rightfully be viewed with caution. Sk 143 is bright enough that follow-up observa- 
tions should be able to significantly reduce the statistical uncertainties without much difficulty. We 
have not pushed the limits of UVES with the present observations. At the same time, measurement 
of A(^Li) in sub-solar metallicity environments are potentially incredibly important. A significant 
abundance of ^Li in low metallicity environments could come from pre-galactic cosmic rays ac- 
celerated via the collapse of large scale structure (producing Li via « + a fusion)^^ or even from 
the effects of non-Standard Model particles in the early universe that affect BBN."' Future obser- 
vations of the gas along this and other SMC or LMC sight lines may be the best way to reliably 
constrain the Li isotopic ratio at sub-solar metallicities. 

6 Notes on Chemical Evolution Models 

We compare our results in Figures [2] and [3] with those of recent chemical evolution models of 
Prantzos,*^ although there are other models in the literatureJ ^^*^^* ^ These models account for post- 
BBN contributions to the total ^Li abundance from massive stars through neutrino nucleosynthesis 
in core collapse supernovaeP' and from low mass stars as red giant or asymptotic giant branch stars 
or novacf^The ^Li yields from these sources are poorly known. Both ^Li and ^Li are synthesized 
directly in the ISM through cosmic ray nucleosynthesis. 

The Prantzos models have worked around the uncertainties in the stellar yields by assuming 
the low-mass star contributions are sufficient to reproduce the solar system (meteoritic) abun- 
dances. The relative contributions of the various sources as a function of metallicity follow the 
relative importance of core collapse (high mass) and type la (low mass) supernovae (SNe). The 
cosmic ray model adopted follows from an analysis of the the energy input from core collapse SNe 
of massive stars that have suffered mass loss through stellar winds.l^ In assessing the yields, no 
credence is given to the Li abundances of stars in the Milky Way. The only data points matched by 
the models are the primordial and meteoritic abundances of Li. 

While the uncertainties in the yields are significant, these models are illustrative of the shape 
of the post-BBN evolution of cosmic Li. However, the primordial and solar system abundances 
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adopted by Prantzos are different than those adopted in our work. Thus, for comparison with 
our measurements, we have had to make adjustments to the model results. We have done this by 
adopting the Cyburt et al.^ primordial abundance, assuming the contributions from core-collapse 
SNe and cosmic rays are as given in Prantzos' work, and independently scaling his adopted low- 
mass star production rates to match our adopted solar system abundance.!^ 



29. Sofia, U.J., Gordon, K.D., Clayton, G.C., Missek, K., Wolff, M.J., Cox, N.L.J., & Ehrenfre- 
und, P. Probing the Dust Responsible for Small Magellanic Cloud Extinction. Astrophys. J., 
636, 753-764 (2006) 

30. Lehner, N., & Howk, J. C. Highly ionized plasma in the Large Magellanic Cloud: evidence 
for outflows and a possible galactic wind. Mon. Not. R. Astron. Soc, 377, 687-704 (2007) 

31. Lehner, N., Howk, J. C, Keenan, F. P., & Smoker, J. V. Metallicity and Physical Conditions in 
the MageUanic Bridge. Astrophys. J., 678, 219-233 (2008) 

32. Gordon, K. D., Clayton, G. C, Missek, K. A., Landok, A. U., & Wolff, M. J. A Quantka- 
tive Comparison of the Small Magellanic Cloud, Large Magellanic Cloud, and Milky Way 
Ultraviolet to Near-Infrared Extinction Curves. Astrophys. J., 594, 279-293 (2003) 

33. Fitzpatrick, E. L. Ultraviolet interstellar absorption toward stars in the Small Magellanic 
Cloud. Ill - The structure and kinematics of the Small Magellanic Cloud. Astrophys. J. SuppL, 
59, 77-93 (1985) 

34. Bouchet, P., Lequeux, J., Maurice, E., Prevot, L., & Prevot-Burnichon, M. L. The visible and 
infrared extinction law and the gas-to-dust ratio in the Small Magellanic Cloud. Astron. & 
Astrophys., 149, 330-336 (1985) 

35. Stanimirovic, S., Staveley-Smith, L., Dickey, J. M., Sault, R. J., & Snowden, S. L. The large- 
scale HI structure of the SmaU MageUanic Cloud. Mon. Not. R. Astron. Soc, 302, 417-436 
(1999) 



19 



36. Stanimirovic, S., Staveley-Smith, L., & Jones, P. A. A New Look at the Kinematics of Neutral 
Hydrogen in the Small Magellanic Cloud. Astrophys. J., 604, 176-186 (2004) 

37. Welty, D. E., & Crowther, P. A. Interstellar Till in the Milky Way and Magellanic Clouds. 
Mon. Not. R. Astron. Soc, 404, 1321-1348 (2010) 

38. Morton, D. C. Atomic Data for Resonance Absorption Lines. III. Wavelengths Longward of 
the Lyman Limit for the Elements Hydrogen to Gallium. Astrophys. J. SuppL, 149, 205-238 
(2003) 

39. Webb, J. K. Quasar Absorption Lines. Ph.D. Thesis, (1987) 

40. Murphy, M. T. Probing variations in the fundamental constants with quasar absorption lines. 
Ph.D. Thesis, (2002) 

41. Savage, B. D., & Sembach, K. R. The analysis of apparent optical depth profiles for interstellar 
absorption lines. Astrophys. J., 379, 245-259 (1991) 

42. Howk, J. C, Sembach, K. R., & Savage, B. D. Ionized Gas in the First 10 Kiloparsecs of the 
Interstellar Galactic Halo. Astrophys. J., 586, 249-267 (2003) 

43. Trundle, C, Dufton, P. L., Hunter, I., Evans, C. J., Lennon, D. J., Smartt, S. J., & Ryans, 
R. S. 1. The VLT-FLAMES survey of massive stars: evolution of surface N abundances and 
effective temperature scales in the Galaxy and Magellanic Clouds. Astron. & Astrophys., 471, 
625-643 (2007) 

44. Bouret, J.-C, Lanz, T, Hillier, D. J., Heap, S. R., Hubeny, I., Lennon, D. J., Smith, L. J., & 
Evans, C. J. Quantitative Spectroscopy of O Stars at Low Metallicity: O Dwarfs in NGC 346. 
Astrophys. J., 595, 1182-1205 (2003) 

45. Kom, A. J., Becker, S. R., Gummersbach, C. A., & Wolf, B. Chemical abundances from 
Magellanic cloud B stars. Astron. & Astrophys., 353, 655-665 (2000) 

46. Luck, R. E., Moffett, T. J., Barnes, T. G., Ill, & Gieren, W. P Magellanic Cloud Cepheids - 
Abundances. Astron. J., 115, 605-634 (1998) 

20 



47. Jenkins, E. B. A Unified Representation of Gas-Phase Element Depletions in the Interstellar 
Medium. Astrophys. J., 700, 1299-1348 (2009) 

48. Cartledge, S. I. B., Lauroesch, J. T., Meyer, D. M., & Sofia, U. J. The Homogeneity of Inter- 
stellar Elemental Abundances in the Galactic Disk. Astrophys. J., 641, 327-346 (2006) 

49. Savage, B. D., & Sembach, K. R. Interstellar Abundances from Absorption-Line Observations 
with the Hubble Space Telescope. Ann. Rev. Astron. Astrophys., 34, 279-330 (1996) 

50. Pequignot, D., & Aldrovandi, S. M. V. The ionization balance in H I regions. Astron. & 
Astrophys., 161, 169-176 (1986) 

51. White, R. E. Interstellar lithium - Differential depletion in diffuse clouds. Astrophys. J., 307, 
777-786(1986) 

52. Samland, M. Modeling the Evolution of Disk Galaxies. 11. Yields of Massive Stars. Astrophys. 
J., 496, 155-(1998) 

53. Hunter, L, et al. The VLT-FLAMES survey of massive stars: constraints on stellar evolution 
from the chemical compositions of rapidly rotating Galactic and Magellanic Cloud B-type 
stars. Astron. & Astrophys., 496, 841-853 (2009) 

54. Dufton, P. L., Ryans, R. S. L, Trundle, C, Lennon, D. J., Hubeny, 1., Lanz, T, & AUende 
Prieto, C. B-type supergiants in the SMC: Chemical compositions and comparison of static 
and unified models. Astron. & Astrophys., 434, 1 125-1 137 (2005) 

55. Venn, K. A. A-Type Supergiant Abundances in the Small Magellanic Cloud: Probes of Evo- 
lution. Astrophys. J., 518, 405-421 (1999) 

56. Peimbert, M., Peimbert, A., & Ruiz, M. T. The Chemical Composition of the Small Magellanic 
Cloud H 11 Region NGC 346 and the Primordial Helium Abundance. Astrophys. J., 541, 688- 
700 (2000) 

57. Zhang, H. W., Gehren, T, Butler, K., Shi, J. R., & Zhao, G. Potassium abundances in nearby 
metal-poor stars. Astron. & Astrophys., 457, 645-650 (2006) 

21 



58. Takeda, Y., Zhao, G., Chen, Y.-Q., Qiu, H.-M., & Takada-Hidai, M. On the Abundance of 
Potassium in Metal-Poor Stars. Pub. Astron. Soc. Japan, 54, 275-284 (2002) 

59. Tolstoy, E., Hill, V., & Tosi, M. Star-Formation Histories, Abundances, and Kinematics of 
Dwarf Galaxies in the Local Group. Ann. Rev. Astron. Astrophys., 47, 371-425 (2009) 

60. Baade, D., Cristiani, S., Lanz, T., Malaney, R. A., Sahu, K. S., & Vladilo, G. Reduced up- 
per limits on the equivalent width of interstellar Li I 670.8 towards SN 1987 A. Astron. & 
Astrophys., 251, 253-258 (1991) 

61. Baade, D., & Magain, P. Very low upper limits on the strength of interstellar lithium lines 
toward SN 1987 A. Astron. & Astrophys., 194, 237-242 (1988) 

62. Sahu, K. C, Pottasch, S. R., & Sahu, M. Search for primordial lithium in the interstellar 
medium towards SN 1987 A. Astron. & Astrophys., 207, Ll-(1988) 

63. Vidal-Madjar, A., Andreani, P., Cristiani, S., Ferlet, R., Lanz, T., & Vladilo, G. The interstellar 
spectrum toward SN 1987 A. Astron. & Astrophys., 177, L17-L20 (1987) 

64. Sansonetti, C. J., Richou, B., Engleman, R., Jr., & Radziemski, L. J. Measurements of the res- 
onance lines of ^Li and ^Li by Doppler-free frequency-modulation spectroscopy. Phys. Rev. A, 
52, 2682-2688 (1995) 

65. Welty, D. E., Hobbs, L. M., & Kulkami, V. P. A high-resolution survey of interstellar NA I Dl 
lines. Astrophys. J., 436, 152-175 (1994) 

66. Rosman, K. J. R., & Taylor, P. D. P. Isotopic Compositions of the Elements 1997. Journal of 
Physical and Chemical Reference Data, 27, 1275-1287 (1998) 

67. Knauth, D. C, Federman, S.R., Lambert, D.L., & Crane, P. Newly Synthesized Lithium in the 
Interstellar Medium. Nature, 405 656-658 (2000) 

68. Reeves, H., Fowler, W.A., & Hoyle, F. Galactic Cosmic Ray Origin of Li, Be and B in Stars. 
Nature, 226, 727-729 (1970) 



22 



69. Meneguzzi, M., Audouze, J., & Reeves, H. The production of the elements Li, Be, B by 
galactic cosmic rays in space and its relation with stellar observations. Astron. & Astrophys., 
15, 337-359(1971) 

70. Steigman, G., & Walker, T. P. Production of Li, Be, and B in the early Galaxy. Astrophys. J., 
385, L13-L16 (1992) 

71. Indriolo, N., Fields, B. D., & McCall, B. J. The Implications of a High Cosmic-Ray Ionization 
Rate in Diffuse Interstellar Clouds. Astrophys. J., 694, 257-267 (2009) 

72. Kalberla, P. M. W., Burton, W. B., Hartmann, D., Amal, E. M., Bajaja, E., Morras, R., Poppel, 
W. G. L. The Leiden/Argentine/Bonn (LAB) Survey of Galactic HI. Final data release of 
the combined LDS and lAR surveys with improved stray-radiation corrections. Astron. & 
Astrophys., 440, 775-782 (2005) 

73. Casuso, E., & Beckman, J. E. Production of Lithium in the Galactic Disk. Pub. Astron. Soc. 
Pacific, 112, 942-960 (2000) 

74. D'Antona, F., & Matteucci, F. Galactic evolution of lithium. Astron. & Astrophys., 248, 62-71 
(1991) 

75. Mathews, G. J., Alcock, C. R., & Fuller, G. M. Galactic chemical evolution with low and high 
primordial Uthium. Astrophys. J., 349, 449-457 (1990) 

76. Woosley, S. E., Hartmann, D. H., Hoffman, R. D., & Haxton, W. C. The i/-process. Astrophys. 
J., 356, 272-301 (1990) 



23 



